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Transport in one-dimensional symmetric devices can be activated by the combination of thermal 
noise and a bi-harmonic drive. For the study case of an overdamped Brownian particle diffusing on 
a periodic one-dimensional substrate, we distinguish two apparently different bi-harmonic regimes: 
(i) Harmonic mixing, where the two drive frequencies are commensurate and of the order of some 
intrinsic relaxation rate. Earlier predictions based on perturbation expansions seem inadequate to 
interpret our simulation results; (ii) Vibrational mixing, where one harmonic drive component is 
characterized by high frequency but finite amplitude-to-frequency ratio. Its effect on the device 
response to either a static or a low-frequency additional input signal is accurately reproduced by 
rescaling each spatial Fourier component of the substrate potential, separately. Contrary to common 
wisdom, based on the linear response theory, we show that extremely high-frequency modulations 
can indeed influence the response of slowly (or dc) operated devices, with potential applications in 
sensor technology and cellular physiology. Finally, the mixing of two high-frequency beating signal 
is also investigated both numerically and analytically. 



PACS numbers: 05.60.-k, 07.50.Qx, 87.10.+O 

I. INTRODUCTION 

The spectral density of the response of a system at 
thermodynamical equilibrium to a sinusoidal time mod- 
ulation consists of a delta-like spike centered at the forc- 
ing frequency. In linear response theory, the system re- 
sponse to two sufficiently weak sinusoidal modulations 
with different frequency is well reproduced by the linear 
superposition of the system response to each modulation, 
separately. In other words, the spectral contents of the 
system output coincides with that of the input signal. 
A significant exception is represented by the harmonic 
mixing (HM) of two commensurate input frequencies of 
comparable magnitude Q] : The system response then can 
contain harmonics of both drive frequencies and thus, un- 
der certain conditions, even a dc component. 

In this paper we focus on the regime when at least one 
drive component is characterized by high frequency but 
finite amplitude-to-frequency ratio. Such a fast modula- 
tion modifies the internal dynamics of the system, so that 
its response to the other harmonic drive is sensitive to 
both modulating frequencies, no matter what their ratio. 
Such a frequency coupling, termed vibrational mixing, 
does not fall within the framework of the linear response 
theory, as the amplitude of at least one drive component 
must be appreciably large. Right for this reason, how- 
ever, the regime investigated here is consistent with the 
operating conditions of many real devices and is, indeed, 
of general applicability. 

Our study case is reprensented by a Brownian particle 
moving on a one dimensional substrate subjected to an 
external bi-harmonic force F(t) and a zero-mean valued, 
delta-correlated Gaussian noise £(t). Its coordinate x(t) 
obeys the Langevin equation (LE) 

x = -V'(x)+F(t)+S(t), (1) 



where 

F(t) = Ax cos(Qi* + 4> x ) + A 2 cos{n 2 t + 4> 2 ) (2) 
with A%,A 2 > 0, 

(£(*)> =o, mm) = 2DS{t), (3) 

and V(x) is the periodic potential of a substrate with 
period L = 2n. 

In Sec. II we compare the results of extensive numer- 
ical simulations with earlier perturbation predictions for 
the rectification current (x) /2ir induced by HM. We con- 
clude that, in spite of the abundance of numerical re- 
sults, the analytical description of HM available in the 
literature is still incomplete and, to some extent, un- 
satisfactory. We then introduce the vibrational mixing 
regime. A high-frequency perturbation pumps energy 
into the system forcing free particle oscillations of am- 
plitude ipQ = Aij tt 2 comparable with the system length- 
scale. In Sec. Ill we demonstrate both numerically and 
analytically that the particle response to an additional dc 
drive is extremely sensitive to the high-frequency pump 
parameter ip$. In Sec. IV we extend our approach to 
investigate the rectification current in a rocked ratchet 
driven by a bi-harmonic force with high and low fre- 
quency components (vibrational rocked ratchet) . Finally, 
in Sec. V we consider the case of a ratchet driven by two 
high-frequency beating harmonic forces. We show that 
in such a limit a vibrational ratchet can be assimilated 
to a pulsated ratchet, where the modulation frequency of 
the substrate amplitude corresponds to the drive beating 
frequency 



II. HARMONIC MIXING 

We know from the literature of the 1970's [l[ that a 
charged particle confined onto a nonlinear substrate is 
capable of mixing two alternating input electric fields of 
angular frequencies f^i and Q 2 ; its response is expected 
to contain harmonics of Oi and Q 2 . As a result, for com- 
mensurate input frequencies, i.e., mQx = nf^, the time 
dependent particle velocity would contain a dc compo- 
nent, too. Such a phenomenon, termed in the later lit- 
erature harmonic mixing, is a rectification effect induced 
by the asymmetry of the applied force. In view of gen- 
eral perturbation arguments, HM was predicted to be of 
the (n + m)th order in the dynamical parameters of the 
system 

@||. Lately, HM was re-interpreted as a man- 
ifestation of the ratchet phenomenon 0, even if no 
substrate asymmetry is required to generate a HM sig- 
nal. 

More recently, the HM mechanism has been inves- 
tigated numerically as a tool to control the transport 
of interacting particles in artificially engineered quasi- 
onedimensional channels @, Q • An interesting variation 
of this problem has been proposed in the context of soli- 
ton dynamics, where the combination of two ac driving 
forces was proven to rectify the motion of a kink-bearing 
chain owing to the inherent nonlinearity of a travelling 
kink @. 

Let us consider, for simplicity, the overdamped 
stochastic dynamics driven by the bi-harmonic force 



with 



F(t) = A x cos(fiit) + A 2 cos(0 2 i) 



V(x) = d(l — cos a;) 



(4) 



(5) 



and fl 2 = 2f2i. A truncated continued fraction expansion 
IH led to conclude that in the regime of low temperature, 
rf > D, the nonvanishing dc component (x) of the parti- 
cle velocity would scale like 



(i) 
D 



M 

2D 



2D' 



(6) 



Quite surprisingly, this result suggests that for small 
drive amplitudes and high substrate barriers, Ax,A 2 <C 
D <C d, the HM signal is negative and independent of d, 
at variance with the numerical results reported in Fig. 
^ Numerical simulation runs for increasing d values re- 
veal a resonant (x(d)} curve. This is not unexpected as 
for d — > (flattening substrate) the zero-mean force 
with (F(t)) — 0, cannot sustain a non-null drift current, 
whereas for d — > oo (high substrate barriers) the inter- 
well activation mechanism gets exponentially suppressed 
and the relevant drift current drops to zero. [The con- 
flicting sign in Eq. is likely to be due to an erroneous 
definition in Ref. fj.] 

The numerical dependence of (x) on the amplitude of 
F(t) is also more complicated than expected from the 
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FIG. 1: (Color online) Transport via HM in the cosine po- 
tential for (j>i = <f>2, A\ — A2, and (a) Q.2 = 2fii, (b) 
Q.2 = 4f2i: (x) versus d. Simulation parameters: fii = 0.01, 
D = 0.2, and A\ = 0.2 (triangles), Ax — 0.4 (squares), and 
Ax — 1.1 (circles). 



perturbation estimate ©■ In Fig. [21 the time average 
of x(t) is plotted versus A = Ax — A 2 at different drive 
frequencies f2 2 = 2f2i. For low drive amplitudes the HM 
signal (x) grows indeed proportional to A 3 , as suggested 
by the scaling law ©, but only for a sufficiently high 
noise level D. 

Moreover, Fig. [21 illustrates another interesting prop- 
erty of rectification by HM: at relatively high ac frequen- 
cies (non-adiabatic regime), the curves (x(A)) develop 
regular oscillations for A > 1 with period and ampli- 
tude roughly proportional to f2i. The details of such a 
non-adiabatic mechanism are explained in Ref. 0: On 
setting A at increasingly high values above the depinning 
threshold of V(x), max{| V'(x)|} = 1, the number of sub- 
strate cells the driven particle crosses during one half- 
cycle, increases by one unit, first to the right and then 
to the left, thus causing one full (x) oscillation at regular 
A increments, AA, proportional to fij. Of course, in the 
adiabatic limit, fii — * 0, these oscillations tend to dis- 
appear with AA. Moreover, shortening the drive period 
or lowering the noise level for A > 1 enhances the above 
modulation effect [jj. Finally, on further increasing A 
the cancellation of the right and the left drifts becomes 
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FIG. 2: (Color online) Transport via HM in the cosine po- 
tential (0 for (f>i = 02, Ai — A2, and (a) Q2 = 2fh, 
(b) r?2 = 4fii: (£} versus Ai. Simulation parameters: (a) 
squares: fli = 0.4, D = 0.2; empty circles: fii = 0.8, D — 0.2; 
triangles: fii = 0.01, £> = 0.2; solid circles: fii = 0.05, 
D = 0.4. Inset: squares: fii = 0.4, D = 0.2; empty circles: 
fii = 0.8, D = 0.2; triangles: fii = 0.4, D = 0.1; (b) squares: 
fii = 0.4, D = 0.2; circles: fii = 0.1, £> = 0.2; triangles: 
fii = 0.01, D = 0.2. In both panels d = 1. 
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FIG. 3: (Color online) Transport via HM in the cosine po- 
tential © for 0i = 02, Ai — A2, and (a) Q.2 = 2fii, (b) 
fi2 = 4fii: (x) versus fii. Simulation parameters: solid sym- 
bols: D = 0.2; empty symbols: D — 0.4; triangles: A\ — 0.4; 
squares: Ai = 0.6; circles: Ai = 1.1; in both panels d = 1. 



more and more efficient; as a result the envelope of the 
(A) oscillations in Fig. [2] decays seemingly inversely pro- 
portional to \f~A. 

An independent perturbation approach Q led to the 
following scaling law for the rectification velocity of a 
Brownian particle JU in a cosine potential 10 subject to 
the harmonic force (0} with VI2 = 2£li 



A, 



<*Wd" 
Qi Id/ V2O1 



A 2 



— — • (7) 



This prediction, that applies under the conditions d <C 
ill *C D, reproduces at least qualitatively both the d — > 
branches of Fig. ^and the f2i — > 00 tails of the curves 
(i(fii)) in Fig. |3J We leave the task of a quantitative 
assessment of Eq. (JTJl for future simulation work. Here we 
limit ourselves to remarking that for large commensurate 
drive frequencies, i.e., £!i = mf7 and ^2 = nQ with 
f2o — ► 00, the HM signal drops sharply to zero. 

We now address two totally different drive regimes that 



cannot be assimilated to the HM phenomenon introduced 
above. Having in mind the drive force Q, we consider 
the limits: 

(i) : f^2 - * 00 with fij and A2/Q2 constant; 

(ii) : Q+ — * 00 with f2_, A1/Q1 and A2/Q2 constant. 

Here we have introduced the shorthand notation Q± = 
i (r2i ± f^)- In (i) the frequency and the amplitude 
of one harmonic component of F(t) are taken both 
large, but with constant amplitude-to-frequency ratio; in 
(ii) the frequency and the amplitude of both harmonic 
components are taken large with constant amplitude-to- 
frequency ratios and, additionally, with constant beating 
frequency 20_ . The response of system Q to an ex- 
ternal drive iQJ in regime (i) or (ii) can be well repro- 
duced within Bleckman's perturbation scheme, or vibra- 
tional mechanics |10|. outlined in the forthcoming Sec. 
IIIII Correspondingly, the ensuing frequency mixing ef- 
fect discussed below will be termed vibrational mixing. 
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III. THE VIBRATIONAL MECHANICS 
SCHEME 

Let us consider the overdamped Brownian particle i ffl l 
with F(t) and defined in Eqs. © and (3J|, respec- 
tively. The periodic substrate potential V{x) has period 
L = 2ir and general form 

oo oo 

V(x) — a„ cos(nx) + 6„ sin(nx), (8) 

n— 1 n— 1 

for an appropriate choice of the Fourier coefficients {a n } 
and {&„}. 

Let us consider for simplicity the regime (i) of Sec. [HJ 
that is we assume that one component of F(t) is slow and 
the other one is fast, say, 12 1 < fi 2 ; then, following the 
approach of Refs. [IlL , we can separate 

— > x(t) + ip(t), (9) 

where, in shorthand notation, from now on x{t) rep- 
resents a slowly time-modulated stochastic process and 
ip(t) is the particle free spatial oscillation 

ip(t) = Vosm(n 2 i + ^2) (10) 

with amplitude "0o = ^W^2- On averaging out ?/;(£) over 
time, the LE for the slow reduced spatial variable x{t) 
can be written as 

x = -V\x) + A x cos(Oii + 4>i) + £(t), (11) 

where 

oo oo 

V{x) = a n Jo(nipo) cos(nx) + b n Jo(mpo) sin(nx). 

n— 1 n— 1 

(12) 

Here, we made use of the identities (sm[nip(t))}t = and 
(cos[ntp(t)])t — Jo(mp ), with Jq(x) denoting the Bessel 
function of 0-order [R| - see also inset of Fig. 0] - and 
((...))( representing the time average of the argument 

This is an instance of the adiabatic elimination of a 
fast oscillating observable ?/>(£), with constant am- 
plitude ipQ. As a result, the slow observable x(t) diffuses 
on an effective, or renormalized potential V(x) driven by 
the slow harmonic in Eq. J5J, alone. We remark that 
V{x) depends on the ratio ipo = A2/SI.2, the amplitude 
of its n-th Fourier component oscillating like \J (nipo)\- 
The adiabatic separation @ for £!i <gc is tenable as 
long as the fast oscillation amplitudes are clearly dis- 
tinguishable with respect to the corresponding Brown- 
ian diffusion 0, that is (ijj(t) 2 )t — \i>l > 2Dt 2 with 
<2 = 27r/S72 or, equivalently, 

In the limit r2 2 — > 00 at constant A2/O.2, the approximate 
LE i|ll[l is expected to be very accurate, regardless of the 
value of D. 



We discuss now a simple application of the vibrational 
mechanics scheme in the presence of a dc drive, i.e., fii = 
and 4>i = 0. The simplest choice for the substrate 
potential is 

V(x) = cosx, (14) 

corresponding to setting a\ — 1 and all the remaining 
Fourier coefficients a n ,b n to zero. The reduced prob- 
lem <|llf) - (ll2f) describes the Brownian diffusion in a wash- 
board potential with variable tilt A\ |15j . 

The observable that best quantifies the response of 
such a system to the dc input A\ is the mobility fi = 
(x) I A\ . In Fig. 0] we compare the simulation data for 
the full dynamics QJ-JSJ against the analytic predictions 
for the static limit of the LE fTT |l - l(T2^ (i.e. when Qi = 0, 
4>i = 0) at increasing ratios A2/Q2 of the ac component 
of F(t). The solid curves on display have been obtained 
by computing the analytic expression (11.51) of Ref. [l^ 
for fi. The agreement between simulation and theory 
is surprisingly close even for noise intensities above our 
threshold of confidence l|13|) . 




FIG. 4: (Color online) Mobility versus A2/Q2 in the dc case, 
fii = and 4>i = 0, for different D. The simulation data 
(dots) have been obtained by integrating numerically the 
LE with V(x) given in Eq. I14H and parameter values: 
Ai = 0.23 and Q2 = 0.1. The solid curves represent the 
corresponding analytic prediction (11.51) of Ref. |l5j| for the 
reduced LE <f 1 ip . Inset: amplitude | Jo(j42/rh)| of V(x) (solid 
curve) compared with the static force Ai. 

The dependence of the system mobility on the high- 
frequency input signal is remarkable: 

(1) The curves ^(^42/^2) exhibit apparent oscillations 
with maxima at fj, = 1. These peaks clearly correspond to 
values of tpo — A2/Q2 for which the amplitude of V{x) = 
Joi^o) cos a; vanishes - see inset of Fig. QJ 

(2) For D = one can apply Eq. (11.54) of Ref. [l|] to 
predict ,i = y/1 - [J (A2/n 2 )/A 1 ]2 for A 1 > J (A 2 /n 2 ), 
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and n = otherwise. As a consequence, the mobility 
curve can develop a number of re-entrant peaks - one for 
each side-peak of | </o(^W^2)| higher than A\ (three in 
Fig. the peak at zero does not count); 

(3) At finite noise intensities, /u, is positive definite for 
any the mobility peaks are located as in the noise- 
less case, but grow less and less sharp as D increases; 

(4) The reduced LE lfH )l -iP J) holds good for A 1 = 0, 
too. This implies that, for an appropriate choice of tpo, a 
high-frequency sinusoidal drive Ai cos(f22i + 02) can neu- 
tralize the effective substrate potential V(x). Explicit nu- 
merical simulations (not shown) substantiate this claim. 
For instance, the time-averaged probability density P(x) 
of the stochastic process , (|2J| and (|14|l flattens out for 
ipo approaching a zero of the Bessel function Jo(V>o)- 

Properties (l)-(4) fully establish the asymptotic regime 
$12 — > 00 at constant A2/SI2 for the dynamics 




2 4 6 8 10 12 14 



A 2 /a 2 

FIG. 5: (Color online) Mobility versus ^2/^2 for the doubly- 
rocked ratchet Q and 1151 1 with Ai — 0.5, Qi = 0.01, cj>i = 
4>2 = 0, and different values of the noise intensity D. All 
simulation data have been obtained for H2 = 10, but the 
black crosses where we set D — 0.12 and Q.2 = 20. Bottom 
inset: simulation data for (1^2/0.2) as in the main panel with 
an additional curve at D — 0.6. Top inset: fi versus D for 
A2 = 0, Ai = 0.5, and Q\ — 0.01; circles: simulation data; 
solid curve: adiabatic formula (11.44) of Ref. [l5| . 



IV. ROCKED VIBRATIONAL RATCHETS 

We consider now a more complicated example that falls 
under the category of the rocked ratchets |T(j . The mo- 
tion of a Brownian particle on an asymmetric substrate 
gets rectified when driven by a time-correlated force, ei- 
ther stochastic or periodic |j| . Let the Fourier coefficients 
of the expansion JHJ all be zero but b\ = — 1 and 62 = —\, 



i.e., 

V(x) = — sin a; — — sin 2a:. (15) 

The corresponding LE (JTJ describes a doubly-rocked 
ratchet 7]. For arbitrary input frequencies f2i, the 
rectified current of the system is known to exhibit marked 
commensuration effects and a complicated dependence 
on the noise intensity and all forcing parameters 0, El • 
We claim here that a well-defined adiabatic limit exists 
for Oi — ► and Q2 — ► 00 with ^2/^2 constant, as sug- 
gested by the separation scheme ©. Following the nota- 
tion of Refs. II ll , we term a rocked ratchet operated 
under such conditions a rocked vibrational ratchet ( VR) . 

The results of our simulation work are summarized in 
Figs. and |SJ To explain the persistent oscillations of 
the curves ^(^/f^), we write down explicitly the renor- 
malized potential i.e. 

V(x) = — Jo(tpo) sinx — - J (2i/> ) sin2x. (16) 

As long as our adiabatic elimination procedure applies, 
the ratchet current j — (x) /2ir vanishes in correspon- 
dence of the zeros of either Bessel function in Eq. i|16fl > 
due to the restored symmetry of the effective substrate. 
On denoting by j n the n-th zero of Jo(x), one predicts 
the following sequence of mobility zeros: 

A 2 1 11 11 , . 

with j a = 2.405, .72 = 5.520, j 3 = 8.654, j 4 = 11.79, 
j 5 = 14.93, etc. [2| . 

As shown in Fig. [S] the sequence (|17|l reproduces 
very closely the zero-crossings of our simulation curves 
for small noise intensities; for D = 0.06 we could locate 
correctly over 20 zeros of the curve ^{A^/Q,^)- In our 
derivation of the effective potential H16|) we cautioned 
that discrepancies may occur for D above the confidence 
threshold l|13|) ; the deviations observed in the bottom in- 
set of Fig. invalidate our approximation scheme only 
for D > 1. The amplitudes of the large ^(^.2/^2) oscil- 
lations decay like (^/f^) - ^ as expected after noticing 
that the modulus of Jq(x) vanishes asymptotically like 
\j2j-Kx for x — > 00 0| . 

Not all zeros of the sequence (|17|l mark an inversion of 
the ratchet current. For instance, for ^2/^2 < \jx the 
current in the effective ratchet potential, (|llfl and l|16ll . 
is certainly positive in the low f2i frequency regime [lfij : 
for i_7i < ^2/^2 < ji, the coefficient of sin 2a; changes 
sign and so does the ratchet polarity (and current); on 
further increasing ^2/^2 larger than ji, the sign of both 
Fourier coefficients Ijl6ll get reversed with respect to l|15ll : 
this is equivalent to turning V(x) upside-down (beside 
slightly re-modulating its profile), so that the polarity 
of V(x) stays negative. Following this line of reasoning 
one predicts double zeros (i.e. no current inversions) at 

■00 = jl,h,j3,ji, ■■■■ 



In the low frequency regime, fli < 1, the reduced 
ratchet dynamics, (fTT|> and can be treated adiabat- 
ically. Its mobility can be com put ed analytically by time 
averaging Eq. (11.44) of Ref. [15j over one forcing cycle 
t\ = 27r/f2i . In Fig. [HI the analytic curves for ^(^2/^2) 
fit very closely our simulation data (dots of the same 
color) at low noise, no matter what the amplitude A\ of 
the slow harmonic in In the bottom inset of Fig. 
deviations from the low frequency curve become visible 
for fix > 0.1: this does not imply that the projection 
scheme leading to the reduced LE tU) - lO fails on in- 
creasing ill with Qi -C ^2) but rather that the adiabatic 
treatment of the reduced LE becomes untenable. This 
conclusion is corroborated by the fact that the mobility 
zeros (and signs) of the curves both in the main panel and 
in the bottom inset of Fig. are independent of either 
parameters A\ and f^i of the low-frequency component. 

Figure [5] illustrates another important VR property. 
In the presence of the high-frequency harmonic, alone, 
A\ = and O2 3> 1, the simulated net current is vanish- 
ingly small (triangles in the main panel). In the absence 
of fast oscillations, A2 = 0, instead, the curve /x(0) versus 
Ai is well reproduced by the adiabatic limit fl± <C 1 0] 
(Fig. H3 top inset). On comparison, one notices that, for 
relatively small Ai, the amplitude of the /^(.^/l^) oscil- 
lations can grow notably larger than the corresponding 
/i(0). This means that energy pumped into the system 
at too high frequency gets dissipated into the heat bath, 
if the system is operated at equilibrium; vice versa the 
nonlinear nature of the system induces a cooperative cou- 
pling between high-frequency disturbances and optimal 
drives, thus enhancing the system response beyond the 
expectations of the linear response theory. 




FIG. 6: (Color online) Mobility versus A2/SI2 for the doubly- 
rocked ratchet Q and with D = 0.12, fii = 0.01, fi 2 = 
10, 4>i = <^2 = 0, and different values of A\. Top inset: fi 
versus Ai for A2 = 0, D = 0.12, and fii — 0.01; circles: 
simulation; solid curves: adiabatic approximation (11.44) of 
Ref. |l5|. Bottom inset: fi versus A2/Q2 for the for the 
doubly-rocked ratchet Q and 115H with A\ — 0.5, D = 0.12, 
Q.2 — 10, 4>i = <j>2 = 0, and different fii. 



so that 



and 



4'1 + "02 



Ai + A 2 

st7 



iln-fazi -^ T (A 1 +A 2 ), 



(20) 



(21) 



PULSATED VIBRATIONAL RATCHETS 



with ipi — tp2 negligible with respect to ipi + ip2 ■ 

On applying the vibrational mechanics scheme of Sec. 
PT1 the effective LE for the reduced spatial variable x(t) 
now reads 



We study now the process 10-© in the regime (ii), 
i.e., for Oi , O2 — > 00 with |02 — ipi = A±/Qi, and 
tp2 = ^2/^2 constant. In order to simplify the algebraic 
passages below, we set <fii = 4>2 > like in the simulations of 
Figs. 17151 

Simple trigonometric manipulations lead to the follow- 
ing expression for the driven free-particle oscillations 



with 



= (ipi + ^2) sin(f2 + i) cos(tt-t) 
+■ (ipi — 1P2) cos(r2+i) sin(f7_t) 

Q + {A 1 ±A 2 )-Q_(A 1 tA 2 ) 



n 2 + 



n 2 



(18) 



(19) 



The parameter range relevant to the discussion of our 
simulation results is 



M +A 2 



A1-A2 



where 

00 



n=l 



x = -V (x,t)+£Qb), 



A1 + A2 
n — — cos(s2_i) 

i £_(_ 



(22) 

[a n cos(nx)+b n sin( 
(23) 



Here we made use of the fact that |f2_| <C f2+, so that 
the time average (cos[nip(t)]) t could be taken over one 
fast oscillation cycle T + = 27r/f2 + , while the slow ampli- 
tude modulation of period T_ = 27r/|57_| was handled as 
an adiabatic perturbation. 

The effect of the bi-harmonic drive F(t) on the sub- 
strate potential is twofold: 

(1) The overall amplitude of the periodic function V(x) 
is modulated in time, which is equivalent to periodi- 
cally modulating the noise intensity D. This condition 
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FIG. 7: (Color online) Transport in a VR Q and GSJ with 
D — 0.06, fi_ = 0.01, (f>i = (j>2 = 0, and different values of 
Ai = A2: (x) versus 



is reminiscent of the so-called temperature ratchets [Tsj 
discussed at length in Ref. Q; 

(2) Since the argument of the Bessel functions in Eq. 
(|23|l depends on the index n, the Fourier coefficients of 
V{x) are distinctly modulated in time; as a consequence, 
their relative weights change in time and so does the pro- 
file ofV{x). 

Separating (sure, rather arbitrarily!) the two time de- 
pendencies (1) and (2) of V{x) helps us explain the sim- 
ulation results of Figs. 17191 obtained for the asymmetric 
potential l|15|) . 
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FIG. 8: (Color online) Transport in a VR JTJ and JT^J with 
Ai = Ai = 0.5, tt~ = 0.01, (f>i = (f>2 — 0, and different values 
of D: (x) versus O+. 

Figures [7| and El illustrate the most interesting feature 
of this class of VR, namely the inversion from positive 



to negative current that takes place for + 3> |0_|. We 
recall that the positive orientation of (x) corresponds to 
the normal polarity of a ratchet (|15|l slowly rocked by a 
harmonic drive, whereas the negative orientation is to be 
expected when the same device is being operated in the 
pulsated mode (like for an adiabatic temperature ratchet 
|18( ). This is why we term regime (ii) of the process 
©-© pulsated VR. 

Note that the pulsated VR of Figs. [7|and|S]are oper- 
ated for values of 0_ and + such that the HM mecha- 
nism plays no significant role, being the HM spikes con- 
fined to values of + much closer to |0_| (as shown in 
Fig. [§). 

The current inversions in the adiabatic regime + 3> 
1 0-| can be explained by looking at the time-dependent 
effective potential 



V(x,t) 



Jo 
1 . 



2A 



• cos(0_i) 



AA 

o~ 



cos(f2_£) 



sin(x) (24) 
sin(2x), 



associated with the substrate potential V(x) of Eq. 115fl . 
Here we set 4>i = (j>2 and Ai = A 2 = A to make con- 
tact with the simulation conditions of Figs. [7| and [S] 
For |4A/n+| < 31, i.e., for fl+ > with Q* + = 4A/j u 
both Fourier coefficients in Eq. (|24|l retain their (nega- 
tive) sign at any time t; V(x, t) does not change polarity 
and the overall effect of the adiabatic modulation with 
period TL amounts to pulsating the amplitude of the ef- 
fective potential (or, equivalently, the noise level [l8|) 
with the same period. As a consequence, (i) is predicted 
to change sign from positive for + < fi^ to negative for 
f2+ > No more current reversals are expected for 

higher f2 + ; as usual in the ratchet phenomenology, (i) 
tends to zero for + — > 00. 

In the actual simulation, see e.g. Figs. 0andEl these 
current inversions seem to take place for f2 + slightly 
smaller than This is due to the fact that at f2+ = 
OI the second coefficient of Eq. (|24|l is negative over 
the entire averaging cycle TL, but for t = mT_ with 
to = 0, 1, 2 ... , where it vanishes. This means that the 
average (x) is still negative at + = 0!j_ and vanishes 
only for lower (but not too lower) + values - in strict 
sense, 0^_ is only an upper bound to the cross-over fre- 
quency. This argument applies as long as condition 
holds; that is not the case of A\ = 0.3 in Fig. [7| and 
D = 1.0 in Fig. 

Finally, we note that the dependence of (x) on the noise 
intensity D exhibits the resonant behavior peculiar to 
most ratchet currents . This happens both for positive 
and negative rectification currents (see Fig. [SJ. 
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FIG. 9: (Color online) From harmonic to vibrational mixing: 
(x) versus fl+ in the bi-harmonically rocked ratchet l|15|l with 
Ai = A2 = 0.5, cj>i = (j>2 = 0, and different values of £1_ = 
0.01 and D. Spikes appear in correspondence with rational 
values of f^/f^i, i.e. for commensurate harmonic drives. 



VI. CONCLUDING REMARKS 

The robustness of the effects simulated here hints at 
the possibility of implementing this concept in the de- 
sign and operation of efficient e.m. wave sensors. As a 
matter of fact, the present investigation has been inspired 
by a typical signal detection problem pj| , namely how to 
reveal a high frequency signal by means of a sensor with 
optimal sensitivity in a relatively low frequency band. 
All cases discussed in the foregoing Sees. IIIIIVI suggest 
a simple recipe: Although the unknown high-frequency 
signal alone cannot be detected, adding a tunable con- 
trol signal with parameters within the device sensitivity 
range causes a nonlinear transfer of energy (information) 
from high to low frequencies, thus enhancing/modulating 
the sensor response to the control signal. By analyzing 
the dependence of the device output on the tunable input 
signal, we can reveal the existence of unknown (and oth- 
erwise undetectable) high-frequency signals. Note that 
in the simulations of Figs. |S] and [5] the input frequencies 



f^i and Q2 differ by four orders of magnitude or more, 
whereas the response function fi can be sensitive to forced 
output oscillation amplitudes ipo as small as the device 
substrate unit length L. In the simulations of Figs. [7|and 
|H1 instead, the ratchet currents are controlled by beating 
frequencies | O2 — two (or more) orders of magnitude 
smaller than the carrier frequency ^(^1 + ^2)- 

As a further application we suggest that the mecha- 
nism of frequency coupling studied here can impact our 
assessment of the health hazards associated with electro- 
pollution |20j| . While high-frequency (non- ionizing) e.m. 
radiation is likely to be harmless at the small length 
scales of sensitive biomolecules, like the DNA helix [2lj . 
still it can affect physiological processes at the cell level. 
Not only high-frequency e.m. waves heat up the biolog- 
ical tissues, but in view the present report, such radia- 
tion can interfere with the much lower-frequency electro- 
chemical control signals that regulate the ion transport 
across cell membranes or the information transfer and 
processing through neuron networks and sensory nerves. 
In conclusion, our results corroborate the recent shift of 
the biomedical research focus from structural (and irre- 
versible) to functional (and possibly reversible) biological 
damages caused by electro-pollution. 

Finally, we mention another potential extension of the 
vibrational mechanics scheme of Sec. IHII In a forthcom- 
ing paper we study the inertial effects of the LE 



x = -yx-V'(x)+F{t)+£{t) 



(25) 



with F(t) and £(*) defined in Sec. [I] and D = jkT. Here 
we limit ourselves to anticipating that the vibrational 
mechanics scheme of Sec. [HI still applies upon replacing 
ipo = A2/O2 with 



^0 — ► ^ 



ipo 



7 2 + ^ 



2 ■ 



(26) 



As a consequence, (x) depends on both ipo an d ^2, which 
suggests the design of highly sensitive devices capable 
of separating the components of a multi-species mixture 
according to the different particle masses. 
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